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Chaos in Partial Differential Equations 

Yanguang (Charles) Li 

Abstract. This is a survey on Chaos in Partial Differential Equations. First 
we classify soliton equations into three categories: 1. (l+l)-dimensional soliton 
equations, 2. soliton lattices, 3. (l+n)-dimensional soliton equations (n > 2). 
A systematic program has been established by the author and collaborators, 
for proving the existence of chaos in soliton equations under perturbations. 
For each category, we pick a representative to present the results. Then we 
review some initial results on 2D Euler equation. 



1. Introduction 

It is well-known that the theory of chaos in finite-dimensional dynamical sys- 
tems has been well-developed. That includes both discrete maps and systems of 
ordinary differential equations. Such theory has produced important mathematical 
theorems and led to important applications in physics, chemistry, biology, and en- 
gineering etc. |H p(| . On the contrary, the theory of chaos in partial differential 
equations has not been well-developed. On the other hand, the demand for such a 
theory is much more stronger than for finite-dimensional systems. Mathematically, 
studies on infinite-dimensional systems pose much more challenging problems. For 
example, as phase spaces, Banach spaces possess much more structures than Eu- 
clidean spaces. In terms of applications, most of important natural phenomena are 
described by partial differential equations, nonlinear wave equations, Yang-Mills 
equations, and Navier-Stokes equations, to name a few. 

Nonlinear wave equations are the most important class of equations in nat- 
ural sciences. They describe a wide spectrum of phenomena; motion of plasma, 
nonlinear optics (laser), water waves, vortex motion, to name a few. Among these 
nonlinear wave equations, there is a class of equations called soliton equations. This 
class of equations describes a variety of phenomena. In particular, the same soliton 
equation describes several different phenomena. For references, see for example |j| 
W. Mathematical theories on soliton equations have been well developed. Their 
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Cauchy problems are completely solved through inverse scattering transforms. Soli- 
ton equations are integrable Hamiltonian partial differential equations which are the 
natural counterparts of finite-dimensional integrable Hamiltonian systems. 

To set up a systematic study on chaos in partial differential equations, we 
started with the perturbed soliton equations. We classify the perturbed soliton 
equations into three categories: 

1. Perturbed (l+l)-Dimensional Soliton Equations, 

2. Perturbed Soliton Lattices, 

3. Perturbed (1 + n)-Dimensional Soliton Equations (n > 2). 

For each category, we chose a candidate to study. The integrable theories for every 
members in the same category are parallel, and for members in different categories 
are substantially different. The theorem on the existence of chaos for each candidate 
can be parallelly generalized to the rest members of the same category. 

• The candidate for Category 1 is the perturbed cubic focusing nonlinear 



Schrodinger equation |22| [21 



14 



id t q = d 2 q + 2[\q\ 2 - cu 2 ]q + Perturbations, 

under periodic and even boundary conditions q(x + 1) = q(x) and q(—x) = 
q(x), w is a real constant. 

The candidate for Category 2 is the perturbed discrete cubic focusing non- 



linear Schrodinger equation [11 23 



1 . 

«<7n = T-j [Qn+1 - %qn + <?n-lj 

+ \ln\ 2 (qn+i + q n -i) ~ 2uj 2 q n + Perturbations , 

under periodic and even boundary conditions q n +N — q-n and q- n = q n . 
The candidate for Category 3 is the perturbed Davey-Stewartson II equa- 



tions [15 



id t q = [dl - d 2 ]q + [2(\q\ 2 - to 2 ) + u y ]q + Perturbations, 

{d 2 + d 2 }u = -4d y \q\ 2 , 
under periodic and even boundary conditions 

q(t, x + l x ,y) = q(t, x, y) = q(t, x,y + l y ) , 
u(t, x + l x , y) = u(t, x, y) = u(t, x,y + l y ) , 

and 

q(t, -x, y) = q(t, x, y) = q(t, x, -y) , 
u(t, -x, y) = u(t, x, y) = u(t, x, -y) . 

We have established a standard program for proving the existence of chaos in per- 
turbed soliton equations, with the machineries: 

1. Darboux Transformations for Soliton Equations. 

2. Isospectral Theory for Soliton Equations Under Periodic Boundary Condi- 
tion. 

3. Persistence of Invariant Manifolds and Fenichel Fibers. 

4. Melnikov Analysis. 

5. Smale Horseshoes and Symbolic Dynamics Construction of Conley-Moser 
Type. 
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The most important implication of the theory on chaos in partial differential 
equations in theoretical physics will be on the study of turbulence. For that goal, we 
chose the 2D Navier-Stokes equations under periodic boundary conditions to begin 



a dynamical system study on 2D turbulence. Since they possesses Lax pair 17] and 



Darboux transformation [25], the 2D Euler equations are the starting point for an 
analytical study. The high Reynolds number 2D Navier-Stokes equations are viewed 
as a singular perturbation of the 2D Euler equations through the perturbation pa- 
rameter e = l/Re which is the inverse of the Reynolds number. Corresponding 
singular perturbations of nonlinear Schrodinger equation have been studied in ]31| 



p0| [19 [EOj. We have studied the linearized 2D Euler equations and obtained a 



complete spectra theorem [16]. In particular, we have identified unstable eigenval- 
ues. Then we found the approximate representations of the hyperbolic structures 
associated with the unstable eigenvalues through Galerkin truncations [18 



2. Existence of Chaos in Perturbed Soliton Equations 

By existence of chaos, we mean that there exist a Smale horseshoe and the 
Bernoulli shift dynamics for certain Poincare map. For lower dimensional systems, 
there have been a lot of theorems on the existence of chaos || J2(| . For perturbed 
soliton equations under dissipative perturbations, we first establish the existence 
of a Silnikov homoclinic orbit. And then we define a Poincare section which is 
transversal to the Silnikov homoclinic orbit, and the Poincare map on the Poincare 
section induced by the flow. Finally we construct the Smale horseshoe for the 
Poincare map. In establishing the existence of the Silnikov homoclinic orbit, we 
need to build a Melnikov analysis through Darboux transformations to generate 
the explicit representation for the unperturbed heteroclinic orbit, the isospectral 
theory for soliton equations to generate the Melnikov vectors, and the persistence 
of invariant manifolds and Fenichel fibers. We also need to utilize the properties of 
the Fenichel fibers to build a second measurement inside a slow manifold, together 
with normal form techniques. The Melnikov measurement and the second mea- 
surement together lead to the existence of the Silnikov homoclinic orbit through 
implicit function arguments. In establishing the existence of Smale horseshoes for 
the Poincare map, we first need to establish a smooth linearization in the neigh- 
borhood of the saddle point (i.e. the asymptotic point of the Silnikov homoclinic 
orbit). Then the dynamics in the neighborhood of the saddle point is governed by 
linear partial differential equations which are explicitly solvable. The global dy- 
namics in the tubular neighborhood of the Silnikov homoclinic orbit away from the 
above neighborhood of the saddle point, can be approximated by linearized flow 
along the Silnikov homoclinic orbit due to finiteness of the passing time. Finally we 
can obtain a semi-explicit representation for the Poincare map. Then we establish 
the existence of fixed points of the Poincare map under certain except-one-point 
conditions. And we study the action of the Poincare map in the neighborhood of 
these fixed points, and verify the Conley-Moser criteria to establish the existence 
of Smale horseshoes and Bernoulli shift dynamics. 

2.1. Existence of Chaos in Perturbed (l+l)-Dimensional Soliton Equa- 
tions. For this category of the perturbed soliton equations, we chose the candidate 
to be the perturbed cubic nonlinear Schrodinger equation. The cubic nonlinear 
Schrodinger equation describes self-focusing phenomena in nonlinear optics, deep 
water surface wave, vortex filament motion etc.. Recently, more and more interests 
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are on perturbed nonlinear Schrodinger equations describing new nonlinear optical 
effects, for example, the works of the Laser Center at Oklahoma State University. 
2.1.1. Dissipative Perturbations. In a series of three papers [E2] 1 21] |14|, we 



proved the existence of chaos in the cubic nonlinear Schrodinger equation under 
dissipative perturbations. We study the following perturbed nonlinear Schrodinger 
equation: 



(2.1) 







M 2 -- 2 


q + ie 



aq + D 2 q + T 



under even periodic boundary conditions 



q(-x)=q(x), q{x + 1) = q(x); 



where i = v— T, q is a complex- valued function of two variables (x, t), (to, a, T) are 
positive constants, e is the positive perturbation parameter, D 2 is a "regularized" 
Laplacian specifically defined by 



D 2 q = 



OO 

E 

3=1 



f3jkj qj cos kjX, 



in which kj = 2nj , qj is the Fourier transform of q, (3j = f3 for j < N, f3j = a*kj 2 
for j > N, [3 and a* are positive constants, and N is a large fixed positive integer. 

Theorem 2.1 (Homoclinic Orbit Theorem). There exists a positive number eo 
such that for any e € (0, eo); there exists a codimension 1 hypersurface E £ in the ex- 
ternal parameter space {tu, a, T, f3, a*}. For any external parameters (u>, a, T, [3, a*) € 
E e , there exists a symmetric pair of homoclinic orbits hk — h)-(t,x) ( k — 1,2) in 
Hi „ ( the Sobolev space H 1 of even and periodic functions ) for the PDE ( 2A ), which 



are asymptotic to a fixed point q e . The symmetry between hi and /i 2 is reflected by 
the relation that hi is a half-period translate of hi, i.e. h 2 {t,x) — h±(t,x + 1/2). 
The hypersurface E e is a perturbation of a known surface [3 — n(u>)a, where k(uj) 
is shown in Figure ^. 



For the complete proof of the theorem, see [22| and [21]. The main argument 
is a combination of a Melnikov analysis and a geometric singuah perturbation 
theory for partial differential equations. The Melnikov function is evaluated on a 
homoclinic orbit of the nonlinear Schrodinger equation, generated through Darboux 
transformations. For more details on this, see the later section on the Darboux 
transformations for the discrete nonlinear Schrodinger equation. 

Theorem 2.2 (Horseshoe Theorem). Un der certain generic assumptions for 
the perturbed nonlinear Schrodinger system (2.1), there exists a compact Cantor 
subset A of E (a Poincare section transversal to the homoclinic orbit), A consists 
of points, and is invariant under P (the Poincare map induced by the flow on T,). 
P restricted to A, is topologically conjugate to the shift automorphism \ on four 
symbols 1,2,-1,-2. That is, there exists a homeomorphism 



4> : W^A, 
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Figure 1. The curve of n — k(oj). 



where W is the topological space of the four symbols, such that the following diagram 
commutes: 



A 



A 



For the complete proof of the theorem, see 14 1. The construction of horseshoes 



is of Conley-Moser type for partial differential equations. 

2.1.2. Singular Perturbations. Recently, singular perturbation, i.e. replacing 
D 2 q by d 2 q, has been studied 30] [19] 1 20 . Consider the singularly perturbed 



nonlinear Schrodinger equation, 

(2.2) iq t = q xx + 2[\q\ 2 - uj 2 }q + is[-aq + (3q xx + T] , 

where q = q{t, x) is a complex- valued function of the two real variables t and x, 
t represents time, and x represents space. q(t, x) is subject to periodic boundary 
condition of period 1, and even constraint, i.e. 

q(t, X + 1) = q{t, x) , q(t, -x) = q(t, x) . 

u> is a positive constant, a > 0, (3 > 0, and T are constants, and e > is the per- 
turbation parameter. The main difficulty introduced by the singular perturbation 
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ed 2 is that it breaks the spectral gap condition of the unperturbed system. There- 
fore, standard invariant manifold results will not apply. Nevertheless, it turns out 
that certain invariant manifold results do hold. The regularity of such invariant 
manifolds at e = is controled by the regularity of e * at e = 0. 



Theorem 2.3 (Li, Jl9[). There exists a e > 0, such that for any e G (0, e ), 
there exists a codimension 1 surface E e in the space of (ui,a,f3, T) G R + x R + x 
R + x R + . where to G (tt,2tt)/S, S is a finite subset. For any external parameters 



on the codimension- one surface, the perturbed nonlinear Schrodinger equation (2.i) 
possesses a symmetric pair of homoclinic orbits hk — hk(t, x) (k = 1,2) in C%° p [0, 1] 
(the space of C°° even and periodic functions on the interval [0,1]) , which is 
asymptotic to a saddle fixed point q £ . The symmetry between hi and hi is reflected 
by the relation that hi is a half-period translate of hi, i.e. h2(t,x) — h\(t, a; + 1/2). 
The hypersurface E e is a perturbation of a known surface [3 = n(u))a, where k(uj) 
is shown in Figure [|. 

2.1.3. Hamiltonian Perturbations. The problem on the existence of chaos in the 
cubic nonlinear Schrodinger equations under Hamiltonian perturbations is largely 
open. The right objects to investigate should be "homoclinic tubes" rather than 
"homoclinic orbits" due to the non-dissipative nature and infinite-dimensionality 
of the perturbed system. Transversal homoclinic tubes are objects of large dimen- 
sional generalization of transversal homoclinic orbits. As Smale's theorem indicates, 
structures in the neighborhood of a transversal homoclinic orbit are rich, structures 
in the neighborhood of a transversal homoclinic tube are even richer. Especially 
in high dimensions, dynamics inside each invariant tubes in the neighborhoods of 
homoclinic tubes are often chaotic too. We call such chaotic dynamics "chaos in 
the smair , and the symbolic dynamics of the invariant tubes ll chaos in the large" . 
Such cascade structures are more important than the structures in a neighborhood 
of a homoclinic orbit, when high or infinite dimensional dynamical systems are 
studied. Symbolic dynamics structures in the neighborhoods of homoclinic tubes 
are more observable than in the neighborhoods of homoclinic orbits in numerical 
and physical experiments. When studying high or infinite dimensional Hamiltonian 
system (for example, the cubic nonlinear Schrodinger equation under Hamiltonian 
perturbations), each invariant tube contains both KAM tori and stochastic layers 
(chaos in the small). Thus, not only dynamics inside each stochastic layer is chaotic, 
all these stochastic layers also move chaotically under Poincare maps. 

There have been a lot of works on the KAM theory of soliton equations under 



Hamiltonian perturbations \29\ Igl [28]. For perturbed nonlinear Schrodinger 



equations, we are interested in the region of the phase space where there exist hyper- 
bolic structures. Thus, the relevant KAM tori are hyperbolic. In finite dimensions, 
the relevant work on such tori is that of Graff Q. In infinite dimensions, the author 
is not aware of such work yet. 



In the paper [13|, the author studied the cubic nonlinear Schrodinger equation 



under Hamiltonian perturbations: 

(2.3) iq t = q xx + 2[\q\ 2 - L0 2 ]q + e[ai + 2a 2 q] , 

under even periodic boundary conditions q(— x) = q(x) and q(x + 1) = q(x); where 
i = q is a complex- valued function of two variable s (t, x) , (u>, ati, 012) are real 

constants, e is the perturbation parameter. The system (|2.3|) can be written in the 
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Hamiltonian form: 

iq t = SH/dq, 

f\\q\ i -2u; 2 \q\ 2 -\q x \ 2 }dx, 
Jo 

H x = [ [a^q + q) + a 2 (q 2 + q 2 )]dx. 
Jo 

Definition 2.4. Denote by W^°' a normally hyperbolic center manifold, by 
\y( cu ) an d \y<, cs ) the center-unstable and center-stable manifolds such that = 
\y( cu ) p| \y( cs )^ anc [ by pt the evolution operator of the partial differential equa- 
tion. Let TL be a submanifold in the intersection between the center-unstable 
and center-stable manifolds W^ cu ' and W^ cs \ such that for any point q € TL, 
distanced* 0),^°)} 0, as |i| — > oo. We call TL a transversal homoclinic tube 
asymptotic to W^ c > under the flow F l if the intersection between W^ cu ^ and 
is transversal at TL. Let E be an Poincare section which intersects TL transversally, 
and P is the Poincare map induced by the flow F*; then HnS is called a transversal 
homoclinic tube under the Poincare map P. 



where H = H + eH\ , 

H a = 



Theorem 2.5 (Homoclinic Tube Theorem). There exist a positive constant Bq > 
and a region £ for (a\,a.2,u)), such that for any e € (— Eq, £q) and any (oti, a%, ui) € 
£ , there exists a codimension 2 transversal homoclinic tube asymptotic to a codi- 
mension 2 center manifold W^- C K 



For a complete proof of this theorem, see 
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2.2. Chaos in Perturbed Soliton Lattices. For this category, we chose the 
candidate to be the perturbed cubic nonlinear Schrodinger lattice. 

2.2.1. Dissipative Perturbations. In a series of three papers [11] 23] [24], we 
proved the existence of chaos in the discrete cubic nonlinear Schrodinger equation 
under a concrete dissipative perturbation. 

We study the perturbed discrete cubic nonlinear Schrodinger equation 

1 



iq n = 



(2.4) 



h? 



is 



aq n 



2q n + q n -i 



+ \q n \ (q n +i + q n -i) - 2w q„ 



h 2 



(q 



n+l 



2q n 



q„ 



+ r 



under even periodic boundary conditions (qN-n = Qn) and (q n +N = qn) for ar- 
bitrary N; where i = \J — 1, q n 's are complex variables, h — 1/N, (uj, a, (3,T) are 
positive constants, e is the positive perturbation parameter. 
Denote by £jv (N > 7) the external parameter space, 

S w = |(w,a,/3,r) ujE (iVtan^,iVtan^), 

re (0,i),ae (0,a ),/3e (0,A>); 

where ocq and /3q are any fixed positive numbers. 
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Theorem 2.6. For any N (J < N < oo), there exists a positive number Sq, 
such that for any e S (0,eo), there exists a codimension 1 submanifold E e in £jv; 
for any external parameters (to, a, (3,T) on E £ , there exists a homoclinic orbit as- 
ymptotic to a fixed point q e . The submanifold E e is in an 0(e u ) neighborhood 
of the hyperplane (3 = k a, where k — k(lu; N) is shown in Figures and 
v = 1/2 -S , < S < 1/2. 



Remark 2.7. In the cases (3 < N < 6), n is always negative as shown in Figure 
[| Since we require both dissipation parameters a and [3 to be positive, the relation 
(3 = na shows that the existence of homoclinic orbits violates this positivity. For 
N > 7, k can be positive as shown in Figure f§ When N is even and > 7, there is 
in fact a pair of homoclinic orbits asymptotic to a fixed point q £ at the same values 
of the external parameters; since for even N, we have the symmetry: If q n = f(n, t) 
solves ([D|), then q n = f(n+N/2,t) also solves (Q. When N is odd and > 7, the 
study can not guarantee that two homoclinic orbits exist at the same value of the 
external parameters. 

For the complete proof of this theorem, see p3| ] . 

Theorem 2.8 (Horseshoe Theorem). Under certain generic assumptions for 



the perturbed discrete nonlinear Schrodinger system (2.4), there exists a compact 
Cantor subset A of E (a Poincare section transversal to the homoclinic orbit), A 
consists of points, and is invariant under P ( the Poincare map induced by the flow 
onYi). P restricted to A, is topologically conjugate to the shift automorphism \ on 
four symbols 1,2,-1,-2. That is, there exists a homeomorphism 

(f> : W^A, 

where W is the topological space of the four symbols, such that the following diagram 
commutes: 

W -±> A 



A 



For the complete proof of the theorem, see [24]. 

The unperturbed homoclinic orbits for the discrete nonlinear Schrodinger equa- 



tion 
(2.5) 



1 



+ \Qn\ 2 {<ln+i + q n -i) ~ 2w 2 g„ 



- nW 



was constructed through the Darboux transformations which will be presented be- 
low in details. The discrete nonlinear Schrodinger equation is associated with the 
following discrete Zakharov-Shabat system Q: 

(2-6) Vn +X = Ltftpn, 

(2.7) 
where 

B {*) 



7^ 



z ihq n 
ihq n l/z 

1 — z 2 + 2i\h — h 2 q n q n ^i + uj 2 
-izhq n -\ + {l/z)ihq n 



-zihq n + (l/z)ihq n -i 
l/z 2 — 1 + 2i\h + h 2 q n q n -i - lu 
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Figure 2. The curve of n = n(w; N) . 
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Figure 3. The curve of n = k(w; N) . 
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and where z = exp(zA/i). 

Fix a solution q n (t) of the system (2.5), for which the linear operator L n has 
a double point z d of geometric multiplicity 2, which is not on the unit circle. We 
denote two linearly independent solutions (Bloch functions) of the discrete Lax pair 



(2.e 



2.7) at z — z d by (<j>+,<j> n ). Thus, a general solution of the discrete Lax pair 



2~7| ) at (q„(t), z d ) is given by 

(j) n {t]Z d ,c) 



+ cq 



where c is a complex parameter called Backlund parameter. We use 
transformation matrix r n by 



to define a 



r, 



z + (l/z)a n 



-l/z + zd n 



where, 



{z d ) 2 A, 
1 



|0„2| 2 + |^| 2 |^l| 2 
2 i \„d\2\i |2 



Z d A 

\z d \ 4 - 1 
(z d ) 2 A n 
|z d | 4 - 1 
z d z d A n 
1 



\4>- 



+ z u 



I>nl<l>n2, 



\c/>nl\ 2 + \z d \ 2 \<j> n 2 



Then we define Q n and ^ by 

i 

(2-8) = — — 0-n+lln 

h 

and 

(2.9) 9 n (t; z) = r„(z; z d ; 0„)^ n (i; z) 

where solves the discrete Lax pair (2j|2.7) at (q n (t),z). Formulas ( |2.S| ) and ( |2.9D 
are the Backlund-Darboux transformations for the potential and eigenfunctions, 
respectively. We have the following theorem [11]. 

Theorem 2.9 (Backlund-Darboux Transformations). Let q n (t) denote a solu- 
tion of the system ( |g.^ ), for which the linear operator L n has a double point z d of 
geometric multiplicity 2, which is not on the unit circle and which is associated with 
an instability. We d enote tw o linearly independent solutions (Bloch functions) of 
the disc rete L ax p air (2.t.~2/\) at (q n , z d ) by (</>+, </>~). We define Q n (t) and ^ n (t; z) 
by (§[|j and fiUl) . Then 

1. Q n (t) is also a solution of the system (2.1). (The eveness of Q n can be 
guaranteed by choosing the complex Backlund parameter c to lie on an certain 
curve.) 

2. ty n (t;z) solves the discrete Lax pair (2.t.2. r i) at (Q n (t),z). 

3. A(z;Q n ) = A(z;q n ), for all z S C, where A is the Floquet discriminant. 



12 



YANGUANG (CHARLES) LI 



4. Q n (t) is homoclinic to q n (t) in the sense that Q n (t) — > e l6± q n {t), expo- 
nentially as exp(— o~\t\) as t — > ±00. Here 9± are the phase shifts, a is 
a nonvanishing growth rate associated to the double point z d , and explicit 
formulas can be developed for this growth rate and for the phase shifts 6± . 

Next we consider a concrete example. Let 
(2.10) q n = q,Vn; q = aexp{-2i[(a 2 - oj 2 )t] + i-y}, 

where N tan £ < a < iVtan ^ for N > 3, 3 tan § < a < 00 for N = 3. Then Q n 
defined in ( |2.S| ) has the explicit representation: 

G 



(2.11) 
where 



iV,w,7,r,±) = g 



1 



G = 1 + cos 2P - zsin2Ptanhr, 
1 



H n = l± 



costf 



sinP sech rcos2ni9, 



t = 47V 2 ^psini9Vp cos2 # - 1 * 
where r is a real parameter. Furthermore, 

v/ p cos 2 1) — 1 

P = arctan — — ; — - — , 

y/psm v 



As t — > ±00, Q„ — * qe Tl2P . Therefore, Q„ is homoclinic to the circle |<jr n | = a, and 
heteroclinic to points on the circle which are separated in phase of — 4P. 

2.2.2. Hamiltonian Perturbations. In the paper [12|, the author studied the 
discrete nonlinear Schrodinger equation under Hamiltonian perturbations: 



(2.12) iq v 



1 



Tafan+l ~ 2< ?« + 1n-l] + |<?n| (Qn+1 + Qn-l) - 2 W q„ 



+ £<{ [ai(g„ + q n ) + a 2 {q n + q n )]q n + [ai + 2a 2 q n ]j^ hl/5„ 



where i = y— T, (j^s are complex variables, n E Z, (u>, 0.1,0(2) are real constants, 
£ is the perturbation parameter, h is the step size, h = l/N, N> 3 is an integer, 
p n = 1 + h 2 \q n \ 2 , and q n +N = qn , q-n = qn- The system ( 2.12 ) can be written 
in the Hamiltonian form: 

dH 

iqn = Pn 



dq n ' 



where H = H + eH\ , 



1 



N-l 
\ " 



n=0 



Pi 



n=0 



"2(^ + g£)] In Pr, 
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Theorem 2.10 (Homoclinic Tube Theorem). There exist a positive constant 
Eq > and a region £ for (ai,ct2,ui), such that for any e G (— £o,£o) a nd any 
(ai,ot2,u>) S £, there exists a codimension 2 transversal homoclinic tube asymptotic 
to a codimension 2 center manifold W^ c \ 



For a complete proof of this theorem, see 
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2.3. Chaos in Perturbed (1 + rt)-Dimensional Soliton Equations (n > 

2). For this category of the perturbed soliton equations, we chose the candidate to 
be the perturbed Davey-Stewartson II equations. The Davey-Stewartson II equa- 
tions describe nearly one-dimensional water surface wave train Q . There have been 
a lot of studies on the inverse scattering transforms for this set of equations jjj] Q . 
The inverse scattering transforms for (l+n)-dimensional soliton equations (n > 2) 
are substantially different from those for (l+l)-dimensional soliton equations and 
soliton lattices. In fact, the Davey-Stewartson II equations possess finite-time sin- 



gularities [27]. For the perturbed Davey-Stewartson II equations, the theory on 
chaos is largely unfinished. So far, its Melnikov theory has been successfully built. 

Although the inverse spectral theory for the DSII equations is very different 
from those for (l-l-l)-dimensional soliton equations and there is no Floquet spec- 
tral theory, its Backlund-Darboux transformation is as simple as those for (1+1)- 
dimensional soliton equations, e.g. the cubic nonlinear Schrodinger equation. These 
Backlund-Darboux transformations are successfully utilized to construct hetero- 
clinic orbits of Davey-Stewartson II equations through an elegant iteration of the 
transformations. In ]2g] , we successfully built Melnikov vectors for the focusing cu- 
bic nonlinear Schrodinger equation with the gradients of the invariants Fj defined 
through the Floquet discriminants evaluated at critical spectral points. The invari- 
ants Fj's Poisson commute with the Hamiltonian, and their gradients decay expo- 
nentially as time approaches positive and negative infinities - these two properties 
are crucial in deriving and evaluating Melnikov functions. Since there is no Floquet 
discriminant for Davey-Stewartson equations (in contrast to nonlinear Schrodinger 



equations [22]), the Melnikov vectors here are built with the novel idea of replacing 
the gradients of Floquet discriminants by quadratic products of Bloch functions. 
Such Melnikov vectors still maintain the properties of Poisson commuting with the 
gradient of the Hamiltonian and exponential decay as time approaches positive 
and negative infinities. This solves the problem of building Melnikov vectors for 
Davey-Stewartson equations without using the gradients of Floquet discriminant. 
Melnikov functions for perturbed Davey-Stewartson II equations evaluated on the 
above heteroclinic orbits are built. 

2.3.1. Darboux Transformations. First we study the Darboux transformations 
for the Davey-Stewartson II (DSII) equations: 

(2.13) id t q = [dl - d 2 y ]q + p%| 2 - u 2 ) + u v ]q , 



\dl+d 2 ]u = -4dj- 12 



under periodic boundary conditions q(t,x + l x ,y) = q(t,x,y + l y ) = q(t,x,y), where 
q and u are a complex- valued and a real- valued functions of three variables (t, x, y). 
To simplify the study, we may also pose even conditions in both x and y. The DSII 
equations are associated with a Lax pair and a congruent Lax pair. The Lax pair 
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is: 

(2.14) 
(2.15) 

where i\i — (f/'i, ip2) T , and 



Lip = Xip , 



D' 



L 



D- 



A = i 



~d 2 x qd x 
rd x dl 



n (£)+g) 
-(D~r) r 2 



Here we denote by 

(2.16) D+ = ad y + d x , D~ = ad y - d x . 

where r — q, a 2 = — 1, 



1 



1 



n = ^[-U + iV], r 2 = -[U + iV], U = 2(\q\"-L0 2 )+u y . 



The congruent Lax pair is: 

(2.17) 

(2.18) 

where ip = (ipi,ip2) T , and 



Lip = Xip , 



L = 



q 

d- 



A = 



-dl qd x 
rd x dl 



-r 2 -(D-q) 
(D+r) -n 



Let (q, r — q, r±, r%) be a solution to the DSII equation, and let Ao be any value of 
A. Denote by ip — (ipi,ip 2 ) T the eigenfunction solving the Lax pair (2.14, 2.15) at 
(<7, r = q, ri, r 2 ; Ao). Define the matrix operator: 

A + a b 
c A + d 

where A = ad y — A, and a, b, c, d are functions defined as: 



L = 



— [ip2 A 2 i>2 + Pipi Ai ipi] , 



b = -^[ip2 Ai tpi - ipi A 2 1P2] 



— [ipi Ai ip 2 ~ 1P2 A 2 ipi] 



d = — [ip 2 Ai 1P2 + 04>i A 2 ipi] , 
in which Ai = ad y — Ao, A 2 = ad y + Ao, and 

A = - [/?|Vi| 2 + |V 2 | 2 ] • 
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Define a transformation as follows: 

(q,r = f3q, r l ,r 2 ) 



(Q, R, R 2 ) . 



(2.19) 



Q 

R 

Rn 



q-2b, 
(3q-2c, 
r l +2{D+a) 
r 2 - 2{D-d) 



$ = T0; 

where <j> is an eigenfunction solving the Lax pair ( [2.14 , [2.15] ) at (q, r — q, ri,7"2; A), 
D+ and D~ are defined in fl2.16| ), 

Theorem 2.11 (|15|). The transformation (2A^j is a Bdcklund-Darboux trans- 
formation. That is, the functions (Q,R. = Q, R\, R2) defined through the transfor- 
mation ( 2.1{\ ) are also a solution to the Davey-Stewartson II equations. The func- 
tion $ defined through the transformation (2.H ) solves the Lax pair (2.14, 2. It ) at 
(Q, R—Q, R\,R 2 ; A). 

A concrete example with two iterations of the Darboux transformations has 
been worked out in 
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2.3.2. Melnikov Vectors. The DSII equations can be put into the Hamiltonian 
form, 



(2.20) 
where 



iqt 
iq t 



SH/Sq , 
-SH/Sq 



H = 



Jo 



[k«r 



kx\ 2 + \i r 2 - n) \q\ 2 ]dxdy. 



Let ip = {ipi,'4> 2 ) T be an eigenfunction solving the Lax pair ( [2.14 , 2.15), and = 
C0i, "02) T be an eigenfunction solving the corresponding congruent Lax pair ( [2.17 , 
[2.18D ; then 



Lemma 2.12. The inner product of the vector 



U = 



s 



^2^2 



whe 



S = 



1 

1 



(2.2L ) vanishes, 



where 



with the vector field JV ' H given by the right hand side of 
(U , JVH) = . 



(/,<?> 



JO 



{/i5i + /2S2} dxdy. 
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and 



1 

1 



Consider the perturbed DSII equations 

i m nil — ' ' J ~ ' 1 '"' ! ' 

(2.21) 



f id t q = [d 2 x - dl\q + [2(M 2 - c 2 ) + u y ]q + eif , 

\ [dl + d$U = -Ad y \q\ 2 , 

where / is the perturbation which can depend on q and q and their derivatives and 
t, x and y. Let G = (f, f) T . Then the Melnikov function has the expression, 



M 



(U, G) dt 



(2.22) 



= 2 



OO pi 



oo JO JO 



R, 



rft. 



where the integrand is evaluated on an unperturbed heteroclinic orbit obtained 
through the Backlund-Darboux transformations given in Theorem 2.11. A concrete 
example has been worked out in 
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3. Two-Dimensional Euler Equations 

One of the most important implications of chaos theory of partial differential 
equations in theoretical physics will be on the study of turbulence. For that goal, the 
author choose the 2D Navier-Stokes equations under periodic boundary conditions 
to begin a dynamical system study. 



(3.1) 



on 

~dt 

On 

dx 



dv 

dy 



on 

dx 



= 0; 



dy 



under periodic boundary conditions in both x and y directions with period 2ir, 
where Q is vorticity, u and v are respectively velocity components along x and y 
directions, e = 1/Re, and / is the body force. When e — 0, we have the 2D Euler 
equations, 

on dn dfl 

~dt 



(3.2) 



dx dy ' 



du dv 
dx dy 



The relation between vorticity f2 and stream function \I/ is, 

_ dv du 

= — — — = a* , 

dx dy 

where the stream function ^ is defined by, 

ay _ a* 

dy ' dx 
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3.1. Lax Pair and Darboux Transformation. The main breakthrough in 
this project is the discovery of the Lax pair for 2D Euler equation [17]. The 



philosophical significance of the existence of a Lax pair for 2D Euler equation is 
beyond the particular project undertaken here. If one defines integrability of an 
equation by the existence of a Lax pair, then 2D Euler equation is integrable. 
More importantly 2D Navier-Stokes equation at high Reynolds numbers is a near 
integrable system. Such a point of view changes our old ideology on Euler and 
Navier-Stokes equations. 

Starting from Lax pairs, homoclinic structures can be constructed through 
Darboux transformations 
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Indeed, in |25j], the Darboux transformation for 
the Lax pair of 2D Euler equation has been found. Our general program is to first 
identify the figure eight structures of 2D Euler equation, and then study their conse- 
quence in 2D Navier-Stokes equation. The high Reynolds number 2D Navier-Stokes 
equation is viewed as a singular perturbation of the 2D Euler equation through 
the perturbation eA, where e = 1/Re is the inverse of the Reynolds number. 
As mentioned above, singular perturbations have been investigated for nonlinear 
Schrodinger equations. 



We consider the 2D Euler equation, 



(3.3) 



where the bracket { , } is defined as 
{f,9} = (d x f)(d y g)- 



(d v f)(d x g) , and SI = A* 



(3.4) 



Theorem 3.1 (|17|). The Lax pair of the 2D Euler equation (3.5) is given as 

f L(p = \<p , 
\ d t f + Alp = , 

where 

L<p = {tt ) <p}, A<p = 
and X is a complex constant, and ip is a complex- valued function. 



In [25 , A Backlun d-D arboux transformation is found for the above Lax pair. 
Consider the Lax pair (3.4) at A = 0, i.e. 

(3.5) {O,f>} = 0, 

(3.6) d t p + {V,p} = 0, 

where we replaced the notation ip by p. 



Theorem 3.2. Let f = f(t,x,y) be any fixed solution to the system (3.1, 3.t ) 
we define the Gauge transform Gj: 

(3.7) p = G f p=±-\p x -{d x \nf)p] , 

and the transforms of the potentials Q and ^ : 

(3.8) * = * + F, = + AF , 
where F is subject to the constraints 

(3.9) {O,AF} = 0, {n,F} = 0. 
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Then p solves the system \3.!\ , \3.(\ ) at (fi, ^>). Thus ftS.y and ( jg.^ ) form the 
Darboux transformation for the 2D Euler equation ( [<?-4 i and its Lax pair ( 3.5\ 3.t ). 



3.2. Linearized 2D Euler Equations. Under the periodic boundary condi- 
tion and requiring that both u and v have means zero, 



u dxdy 



v dxdy = 0, 



expanding fl into Fourier series, f2 = X)fcez 2 /fo> L ° k e , where uj-k = ^>k , k — 



(&i,fc2), = (^,y), the system (3.2) can be rewritten as the following kinetic 
system, 



(3.10) 

where A(p, q) is given by 
(3.11) 



fc=p+g 



= 2 -\p\ 2 }(PiQ2 - p 2 qi) 



M" 2 -br 2 ] 



Pi <7i 

P2 q-2 



where \q\ 2 = q\ + q\ for q = ((71,(72), similarly for p. To understand the hyperbolic 




(P2,~PlY 



Figure 4. An illustration of the classes S^, and the disk D\ p \. 



structures of the 2D Euler equations, we first investigate the linearized 2D Euler 
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3{A} 





i2\b\ 









-i2\b\ 



3?{A} 



Figure 5. The spectrum of La in case 1. 
3{A} 



• 


i2\b\ 

• 

• 

• 















• 






-t2|&| 



Figure 6. The spectrum of Ca in case 2. 



equations at a stationary solution. Denote {wk\kez 2 /{o} by u). Consider the simple 
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3{A} 



5R{A} 



Figure 7. The quadruple of eigenvalues for the system led by the 
class Sj, labeled by k — (1, 0) T , when p = (1, 1) T . 



fixed point lo*: 
(3.12) 



0, if k ^ p or — p. 



of the 2D Euler equation (3.10), where T is an arbitrary complex constant. The 
linearized two-dimensional Euler equation at lo* is given by, 

(3.13) uj k = A(p, k-p)T u) k -p + A(-p, k+p)T uj k+p . 

Definition 3.3 (Classes). For any k € Z 2 /{0}, we define the class to be 
the subset of Z 2 /{0}: 



k + npe Z 2 /{0} 



n G Z, p is specified in (3.12) 



See Figure [| for an illustration of the classes. According to the classification 
defined in Definition the linearized two-dimensional Euler equation ( |3.13 ) de- 
couples into infinite many invariant subsystems: 



k+np 



(3.14) 



A(p,k+(n-l)p) rw fc+(n _ x)j , 

+ A(-p, 'k+(n + l)p) f uj k+(n+1)p 



Definition 3.4 (The Disk). The disk of radius \p\ in Z 2 / {0}, denoted by D\ p \, 
is defined as 



D H = \keZ 2 /{0} 



\k\ < \ P \ 
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Figure 8. The heteroclinic orbits and unstable manifolds of the 
Galcrkin truncation. 
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The closure of D\ p \, denoted by D\ p \, is denned as 

D M = {keZ 2 /{0} \k\<\ P \ 



Theor em 3 .5 (Unstable Disk Theorem). If S^nL>| p | = 0, then the invariant 
subsystem ( 3. 14\ ) is Liapunov stable for all t G R, in fact, 



E |"W*) " ^ a E |«W°) 



where 



max{-p„} 



min{-p„} 



Vt e i?, 



< a < oo . 



Theorem 3.6. TTie eigenvalues of the linear system (3.14) are of four types: 
real pairs (c,—c), purely imaginary pairs (id, —id), quadruples (±c±id), and zero 
eigenvalues. 

Theorem 3.7 (The Spectral Theorem). 1. //Eg n D\ p \ = 0, then the en- 
tire li spectrum of the linear operator La ( defined by the right-hand side of 
the invariant subsystem) is its continuous spectrum. See Figure 

2. If R D\ p \ 7^ ; then the entire essential ii spectrum of the linear operator 
La is its continuous spectrum. That is, the residual spectrum of La is empty, 
^t{La) = 0- The point spectrum of La is symmetric with respect to both 
real and imaginary axes. See Figure [fj. 

We can calculate the eigenvalues through continued fractions. Let p = (1, 1) T , 
in this case, only one class Sf labeled by k = (1,0) T has no empty intersection 
with Dipi (the other class labeled by k = (0, 1) T gives the complex conjugate of 
the system led by the class labeled by k = (1, 0) T ). For this class, there is no real 
eigenvalue. Numerical calculation through continued fractions gives the eigenvalue: 

A = 0.24822302478255 + i 0.35172076526520 . 

Thus we have a quadruple o f eige nvalues, see Figure [7] for an illustration. Denote 



by L the right hand side of ( 3.13 ), the spectral mapping theorem holds. 



Theorem 3.8 (|10[). 

a(e tL ) = e t ^ L \t^0. 

Moreover, the number of eigenvalues has a sharp upper bound. Let £ denote 
the number of points q £ Z 2 /{0} that belong to the open disk of radius \p\, and 
such that q is not parallel to p. 



Theorem 3.9 (|10|). The number of nonimaginary eigenvalues of L (counting 
the multiplicities) does not exceed 2£. 
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3.3. Approximate Explicit Representations of the Hyperbolic Struc- 
tures of 2D Euler Equations. From Figure [?], we see that the simple fixed point 
given by p = (1,1), has unstable eigenvalues. Our interest is to obtain represen- 
tations of the correponding hyperbolic structures for 2D Euler equations. In 
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through Galerkin truncation, we obtained the approximate explicit representation. 
Figure || shows the heteroclinic orbits and unstable manifolds of the Galerkin trun- 



cation. 



4. Conclusion and Discussion 

We have reported the status of chaos in nonlinear wave equations and of study 
on 2D Euler equations. In particular, we have summarized the most recent results 
on Lax pair and Darboux transformations for 2D Euler equations. 
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